Branching Transition of a Directed Polymer in Random Medium 
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A directed polymer is allowed to branch, with configurations determined by global energy opti- 
mization and disorder. A finite size scaling analysis in 2D shows that, if disorder makes branching 
more and more favorable, a critical transition occurs from the linear scaling regime first studied by 
Huse and Henley [Phys. Rev. Lett. 54, 2708 (1985)] to a fully branched, compact one. At criticality 
clear evidence is obtained that the polymer branches at all scales with dimension d c and roughness 
exponent £ c satisfying (d c — l)/£c — 0-13 ± 0.01, and energy fluctuation exponent uo c — 0.26 ± 0.02, 
in terms of longitudinal distance. 
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Linear directed polymers in random medium (DPRM) 
are a paradigmatic model jl],^] , of central importance for 
fields like equilibrium interfaces in random ferromagnets 
H or fracture ||], in 2D, and flux lines of high-T c super- 
conductors ||, in 3D. The behavior of DPRM at large 
scales is generally controlled by properties at tempera- 
ture T = 0, where the model reduces to an optimiza- 
tion problem JjLS]. Imagine we assign independently a 
random energy St, to each edge b of a square lattice. At 
T = 0, the optimal DPRM is the linear directed (no over- 
hangs) path II(i), which covers a given distance t parallel 
to its direction, e.g. (1, 1), and minimizes En = X^en ^fa- 
Optimization determines a self-affine geometry of the op- 
timal II and a peculiar scaling of En fluctuations, quite 
universal with respect to different forms of the environ- 
mental disorder JlJ . 

So far, much activity on DPRM concentrated on a 
transition occurring in high enough D between a low-T 
regime with rugged free energy landscape, and an high-T 
one, with smooth phase space In the present Letter 
we discuss a novel critical transition, triggered by disor- 
der alone, and occurring in an appropriate generalization 
of DPRM. We release the constraint of linearity and con- 
sider polymers which can form branches and loops. At 
the same time, we assume that the random energy Sj, 
can take both positive and negative values. We look for 
the branched directed polymer LI (BDPRM) spanning a 
certain longitudinal distance t and optimizing En ( de- 
fined as in the linear case) . Even a very little percentage 
of negative energy bonds suffices now to induce branches 
and loops in the optimal II, at least at small scales. In- 
deed, a single negative bond met along an optimal linear 
path can easily become a branch, if it lowers En once in- 
cluded in LL However, as long as the weight of negative 
energies in the f^-distribution remains low, the optimal 
IT is likely to remain linear at large scales. On the other 
hand, if this weight is very large, one expects branches 
and loops to develop wherever possible, leading to a sort 
of fully ramified, compact structure for II. 



For a specific model, we find here that such two limit- 
ing regimes hold on respective sides of a sharp threshold 
concentration of negative energy bonds. Right at thresh- 
old, a critical transition occurs, with II characterized by a 
branching probability between and 1 at all scales. This 
critical regime, with an infinite hierarchy of loops and 
dangling ends, is the first of this kind met in this field 
and demonstrates an unexpected richness of possible so- 
lutions of problems of global optimization with disorder. 
The transition has analogies with theta- or similar points 
in polymer physics. 

A main motivation to study transitions like this comes 
from recent attempts to explain the puzzling universal- 
ity of 2D random Potts ferromagnets [Q. To the pur- 
pose of evaluating its free energy, the interface between 
two phases of a q-state Potts model can be schematized 
as a directed path, which possibly includes bubbles of 
the other q — 2 phases @. At T = these bubbles are 
induced by the presence of antiferromagnetic couplings, 
which naturally proliferate under renormalization. Such 
interfaces with bubbles are ramified polymers without 
dangling ends, and obey global optimization. In the fer- 
romagnetic phase bubbles should not develop at large 
scales, while in the disordered phase they should form 
a compact "froth" structure. The existence and nature 
of peculiar, bubbling critical regimes for these polymers 
can be an important element to understand the above 
universality Another potential field of application 
are models of fracture in disordered materials Q . 

Counting the configurations of a ramified structure, 
even if directed, is a problem of non-polynomial com- 
plexity. An additional difficulty hindering progress in 
these problems comes from the need of sampling ade- 
quately over disorder. Thus, our study here is based on 
a quite extensive computational effort, combined with 
a suitable choice of model and method of scaling analy- 
sis. We consider bond branched polymers directed at 45° 
with respect to the coordinate directions of a square lat- 
tice, with origin in O (Fig. la). No branch can turn back- 
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wards with respect to the polymer direction. A polymer, 
II(t), extends from O up to the straight line, St, at dis- 
tance t from O and orthogonal to the polymer direction 
(Fig. la). For bond energies we choose a distribution 
P(£ b )=p5(£ b - 1) + (l-p)5(£ b + 1) with < p < 1 @. 
Upon decreasing p from 1 , we expect to possibly meet a 
threshold as described above at some p = p c . 



o 




(a) (b) 
FIG. 1. (a) BDPRM configuration without strip restric- 
tion, t is the number of steps of any linear path joining O 
with S t . (b) A portion (At = 21) of optimal BDPRM con- 
fined within a strip (L = 13). Disorder is generated with 
p — 0.925. Strips are periodic in the L direction. 

Information on the optimal path in a given disorder 
configuration is transferred from t to t + 1 by updating 
the set {E(t, t)} of the optimal II energies for given t and 
t = II(t) n St. The recursion for E reads: 

E(T',t + l)=min[E(T,t) + e(a)} (1) 

T,cr 

where t' = H(t + 1) n St+i and e(cr) is the energy associ- 
ated to each admissible set of bonds <t(t,t') connecting 
points in r to points in r'. There are typically several 
such sets for given t and r'. The least energy, E m i n (t) = 
mini?n(t) is then calculated as -E m in(i) = min£(r, t). 

Thus, in order to identify the optimal II, at any given t, 
we have to determine E(t' , t') for all t' < t and r' C St', 
consistently with eq. ([[]). There are 2 t+1 — 1 different r 
intersections corresponding to all non-empty subsets of 
points in St, and altogether 2(i+l) bonds connecting St 
to St+i, hence a transfer algorithm must scan 2 2 * +2 — 1 
cr-configurations to calculate the minimum in (|]). This 
means in total 2 3 ( t+1 ) operations per i-step. Thus, as 
soon as t exceeds a few units, a reasonable statistics over 
disorder can not be collected. Clearly, one cannot simply 
apply the strategy of following the polymer development 
within the triangle in Fig. la up to large t: a different 
method of scaling analysis and some restrictions on al- 
lowed configurations are in order. Thus, we confine II 
within a strip of fixed width, L, and length t (Fig. lb). 
The i-transfer of information on optimal configurations 



within such a strip has a ^-independent computational 
cost determined by L. Our strategy is to extract infor- 
mation on asymptotic properties from strip calculations 
for relatively small L by exploiting finite size scaling ||. 

Since transfer can be made much faster by restrict- 
ing to only one a the possible connections between two 
compatible r and r' in cq.([j]), the following constraint is 
applied to determine such unique a: BDPRM configura- 
tions are assumed to be always such to necessarily imply 
bond connection between neighboring sites in r C St and 
in in r' C St+i- Namely the unique a considered in (Q) is 
that accomplishing all possible nearest neighbor connec- 
tions between the considered r and r' . Such restriction 
offers the advantage of reducing the operations necessary 
to determine E(t', in (P by a factor 2 x 2 L . On the 
other hand, such a represents the contribution to trans- 
fer matrix elements with the highest rate of loop and 
branch formation, leaving also a wide possibility of dan- 
gling ends in the structure, and thus making the model 
fully adequate to display the transition we are interested 
in. Indeed, we also performed extensive calculations for 
an ordered version of this restricted model at T > to 
test the effects of such restriction. On strips with L up to 
10, with methods similar to those in use for directed lat- 
tice animals (DA) ]l0]| , we estimated correlation length 
exponents i/,, = 0.81 ± 0.02 and v ± = 0.48 ± 0.02, to be 
compared with v {l = 0.818 ± 0.001 and v x = 1/2 for DA 
p0| , pl| . This suggests that our model, without disorder, 
belongs to the same universality class as unrestricted DA. 

We identify the transverse BDPRM width, at a given t, 
as the root mean square distance of the points in r from 
the middle of the strip. Since we expect such width, av- 
eraged over random bond configurations, to grow like 
(0 < C 5: 1) on infinite lattice, crossover on our strips 
should be controlled by the ratio t^/L. So, the average 
total number of bonds should scale with L and t as: 

N(L,t)^(N(L,t))=t d n(^j (2) 

where d is a fractal dimension and n is a crossover func- 
tion . The brackets indicate average over disorder. For 
t , L ^> 1, and it <C L we should recover the behavior on 
unrestricted lattice. So lim nix) = const. ^ 0. On the 

other hand, for it L, M must become proportional 
to t as a consequence of statistical invariance under t- 

l-d 

translations. Thus, n(x) ^ x ^oo x 5 . This means that 
Af(L,t) ~ t>i A(L)t, with A(L) - L(^ _1 )/C. In a similar 
way, for the fluctuations of optimal path energies due to 
lattice disorder we expect 



A£(L,t) = ({E min (t) - {E min (t)}) 




For it > L, AS should grow like t ' , as a consequence 
of the central limit theorem |12] . The amplitude of this 
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t 1 ' 2 behavior should scale as B{L) ~ L (u ~ 1/2)/c . On the 
other hand, for < I, e ~ const., as w describes AS 
scaling without strip restriction. Of course, functions like 
n and A, and the various exponents, can depend on p. 

We explored systematically the scalings in eqs.(|J) and 
(||) for different p's. L and t ranged up to 9 and sev- 
eral hundreds, respectively, and data for J\f and AS re- 
sulted from statistics over a minimum of 2 x 10 3 up to a 
maximum of 10 disorder configurations. The most con- 
vincing exponent determinations always occurred when 
analyzing regimes at if- ^> L. E.g., the amplitude A(L,p) 
can always be estimated with good accuracy (relative er- 
ror bars < 10~ 2 ), being determined by averaging over 
large data sets. The trend of the effective exponents 
ln[A(L + l,p)/A(L,p)]/\n[(L + l)/L] can be extrapo- 
lated rather well for all p's (Fig. 2). For p < 0.90 ex- 
trapolation to (d — 1)/C — 1 is always clear. This indi- 
cates a compact, fully branched regime (d = 2, £ = 1). 
On the other hand, for 0.94 < p < 0.98, we extrapolate 
(d — 1)/C ~ |l|. This suggests d = 1, i.e. a linear 
regime. 




FIG. 2. Effective determinations of (d — l)/£ for different 
L and varying p. The insert shows a magnification of the 
transition region, where exponent curves intersect. 

Most interesting, there is a narrow region (p ~ 0.92 
0.93) at which the trends of approach to the above two 
asymptotic values exchange each other, and intersections 
of effective exponent "curves" for different L concentrate 
and stabilize. The corresponding value (d c — 1)/Cc — 
0.12 -r 0.14 is an obvious candidate to represent a dis- 
tinct border regime right at p = p c ~ 0.925. The data in 
Fig. 2 suggest a step at p — p c for L — > oo. The scenario is 
further elucidated by a similar analysis of B amplitudes 
in eq.(|). This time we get clearly (u - 1/2) /C ^ 0.50 
for p <b 0.86 and (u - l/2)/£ ~ -0.27 for 0.94 < p < 
0.98. The value 0.50 confirms a compact, fully branched 



regime, for which we expect £ = 1, and also, rather nat- 
urally, uj = 1 |3- (u; - l/2)/C ~ -0.27 in the d ~ 1 
region is clear evidence of standard, linear DPRM scal- 
ing (w = 1/3 , C = 2/3 (l]J|). Thus, in this regime the 
presence of branchings at small scales does not change 
the universal behavior with respect to the strictly linear 
case. Unfortunately, it is more difficult to guess precisely 
the border-line value of (lu c — l/2)/£ c . This is probably 
due to the circumstance that for p — p c this combination 
of exponents seems to be lower than in the regimes on 
both sides. So, the crossover pattern can not stabilize 
simply and nicely as in the case of Fig. 2. A direct deter- 
mination of u) in the short time regime (t < L = 9) , by 
fitting eq.(||) directly, gives oj c — 0.26±0.02, distinct from 
the DPRM w = 1/3 [UJ. Assuming £ c equal to 1, this 
lu c is consistent with our extrapolation of (u> c — l/2)/£ c 
from B(L,p c ) (uj c = 0.27 ± 0.02). Further insight comes 
from data collapse plots testing A in the form 



A(L 1 p) = L^- 1 >/^a((p~p c )L*) 



(4) 



where <b is a crossover exponent (Fig. 3). Eq.(^) follows 
from (H), once explicitated the dependence on the rele- 
vant parameter (in renormalization group (RG) sense), 
p — p c . Over the range studied, data collapse best for 
p c = 0.927 ± 0.005, (<i c - l)/( c = 0.13 ± 0.01 and 
4> = 1.00±0.02. a(x) — Ae~ ax seems to be an acceptable 
form for the crossover function a. 
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FIG. 3. Collapse plot of the rescaled mass-amplitude data. 

The above analysis shows that, upon increasing the 
fraction of negative bond energies on the lattice, BD- 
PRM undergo a sharp, critical transition from linear 
DPRM, to compact, fully branched regimes. At the tran- 
sition (d c - 1)/Cc = 0.13 ± 0.01, uo c = 0.26 ± 0.02 and 
4> = 1.00±0.02, indicating a fractal regime, with branches 
on all scales, and peculiar energy fluctuations exponent. 
For sure d c , Cc are different from the corresponding expo- 
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nents of DP (d = 1.473 ± 0.001 @, C = 0-6330 ± 0.0008 
0, (d-l)/C = 0.747±0.002), which is a reference model 
for a variety of phenomena, including fracture Jl8| , ^9| . 
Our critical regime appears a peculiar result of optimiza- 
tion and disorder. 

It is natural to ask whether the T = behavior an- 
alyzed above should apply also at T > 0. A numeri- 
cal investigation of T > properties in the random case 
would be computationally too expensive, while higher di- 
mensions are definitely beyond present possibilities. For 
p = 1 and T > our model describes DA with bond 
fugacity exp(— 1/T). Application of Harris' criterion 
|0| to the DA critical point (T c = 0.9708 . . .) suggests 
that disorder is relevant for our model in 2D, because 
z^l [1 + (d— 1)C] < 2 (see exponents quoted above). We 
expect a critical line joining the critical points (1, T c ) and 
(p c ,0) on the (p,T) plane. The point (p c ,0) should be 
attractive from an RG point of view. Indeed, in a coarse- 
grained description, based on a branched path integral, 
for the BDPRM at T > one would expect an elastic 
free energy term -i J (^jf) dt associated with the i-th 
branch. This implies dT/de = (1 — 2QT if t is rescaled by 
a factor 1 + Se and ( is the roughness of a single branch. 
If C > 1/2, which is plausible in our case, irrelevance in 
T follows. So, the whole critical line for p < 1 could be 
controlled by the T = fixed point tested by our scal- 
ing analysis. In higher D, while disorder should remain 
relevant on the basis of Harris' criterion the transi- 
tion of DPRM to smooth high-T free energy landscape 
|| could modify the above phase diagram. 

Summarizing, we gave evidence that properly chosen 
disorder conditions determine a branched critical regime 
in a generalized version of DPRM. In a previous study, 
a similar model without dangling ends on hierarchical 
lattice, displayed a strictly linear (d c = 1) T = Q criti- 
cal regime [Q. So, dangling ends could be essential for 
obtaining d c > 1 in a BDPRM model. Further clarifica- 
tion of this issue may be a crucial test for the scenario 
drawn in ref. fit] , concerning the universality mechanism 
in random ferromagnets with several coexisting phases. 
However, a considerably more substantial computational 
effort would be required to such purpose. A compari- 
son with DP clusters is also appropriate for our critical 
branched polymer. It was shown recently that the DP 
cluster backbone is the structure hosting and shaping (in 
the sense of Hurst's exponent) optimal linear paths for an 
extremal version of the DPRM model . The branched 
fractal cluster identified here, is clearly different from the 
DP one, and originates in a complete and direct way from 
global optimization and disorder. 

We thank Mehran Kardar for useful comments. 
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